Modified two-potential approach to tunneling problems 
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One-body quantum tunneling to continuum is treated via the two-potential approach, dividing the 
tunneling potential into external and internal parts. We show that corrections to this approach can 
be minimized by taking the separation radius inside the interval determined by simple expressions. 
The resulting modified two-potential approach reproduces the resonance's energy and the width, 
both for narrow and wide resonances. We also demonstrate that, without losing its accuracy, the 
two-potential approach can be modified to a form resembling the R-matrix theory, yet without any 
uncertainties related to the choice of the matching radius. 
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I. INTRODUCTION 

The quantum mechanical tunneling through a classi- 
cally forbidden region is an ubiquitous phenomenon in 
physics, vifhich has been extensively studied since the 
early days of quantum mechanics. In 1927, Hund P| 
was the first to point out the possibility of "barrier pen- 
etration" between two discrete states. In the same year, 
Nordheim considered the case of tunneling between 
continuum states. Subsequently, Oppenheimer ^ per- 
formed a calculation of the rate of ionization of the hy- 
drogen atom, and Gamow, 4] Gurncy and Condon Q 
explained alpha decay rates of radioactive nuclei in terms 
of the tunneling effect. 

While the semi-classical treatment of tunneling turned 
out to be very successful in many applications, the nu- 
merical calculation offers very little insight into the phys- 
ical process. In addition, the validity of the standard 
WKB formula is rather restricted. Other methods, al- 
though more accurate, contain various uncertainties. For 
exarnple, the results of the commonly used R-matrix the- 
ory la] are often sensitive to the choice of the matching 
radius 0, , and the theoretical error is difficult to esti- 
mate. 

The treatment of the tunneling problem can be es- 
sentially simplified by reducing it to two separate prob- 
lems: a bound state problem and a non-resonant (scat- 
tering) state problem. This can be done consistently in 
the two-potential approach (TPA) H IH 113 (see also 
Refs. ^2 ^3) J representing the barrier potential as a 
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sum of the "inner" and the "outer" terms, containing 
only bound and only scattering states, respectively. This 
approach not only provides better physical insights than 
many other approximations but it is also simple and ac- 
curate. 

In this paper we propose further developments and a 
modification of the TPA, and present a detailed com- 
parison of this approach with the results of numerical 
calculations based on the Gamow-state (resonant-state) 
formalism. The resulting analytical expressions are easy 
to interpret and they can be straightforwardly extended 
to the non-spherical case. 

The paper is organized as follows. In Sect. II, the TPA 
is briefly described. Section III deals with the quantal 
correction terms to the TPA. The minimization of these 
terms prescribes unambiguously the "window" for the 
separation radius that divides the original barrier poten- 
tial into inner and outer terms. In this case, by consid- 
ering examples of wide and narrow nuclear resonances, 
we demonstrate that the TPA yields results which are 
practically the same as those of the resonant-state cal- 
culation. In Sec. IV, we present a modification of the 
TPA. The resulting expressions resemble those of the R- 
matrix theory, yet without any uncertainties related to 
the matching radius, see Sec. V. Finally, the summary 
of our work is contained in Sec. VI. 



II. TWO-POTENTIAL APPROACH 

Consider a quantum well V{r) with a barrier, which 
contains a quasi-stationary state at the Eres- The coor- 
dinate space can be divided into two regions, the "inner" 
region, < r < i?, and the "outer" region, r>R, where 
R is taken inside the barrier (see Fig. Accordingly, 
one can introduce the two auxiliary potentials: the inner 
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FIG. 1: The "inner" ([/; top) and the "outer" {W; bottom) 
parts of the potential V{r) defined as in Eqs. JTJ and lllj, 
respectively. The separation radius R is chosen well inside 
the barrier. The energy of the metastable state is Eq and 
Vo~V{R). The barrier radius is denoted by R and ri,2 are 
the classical turning points. 



potential 

U{r) = 
and the outer potential 



Vir) 
V{R) 



for r < R 
Vq for r > R 



Vo for r <R 
V{r) for r > R 



(1) 



(2) 



The inner potential contains a bound state, $0(7") 
{E=Eo), representing an eigenstate of the "inner" Hamil- 
tonian Hq = K+U(r), where K — — V^/2m is the kinetic 
energy term {h=l). One can demonstrate ^3 that the 
energy and the width of the quasi-stationary state, asso- 
ciated with the complex-energy poles of the total Green's 
function G{E) = {E ~ K — V)~^ , are obtained from the 
following equation 

E = Eo + ($o|iy|$o> + {<i>o\WG{E)W\^Q) . (3) 



Here W{r) 
given by 



W{r) — Vq, and the Green's function G is 



GiE) = GoiE) 1 + WG{E) 



(4) 



where 



GoiE) 



1 - A 



E + Uo- K -U' 



A=|$o)($o|. (5) 



The resonance energy E'^cs = ^{E) and the width F = 
—2'^{E) of a quasi-stationary state obtained from Eq. Q 
are independent of the choice of the separation radius R. 

Equation (j3Jl can be solved iteratively by using the 
standard Born series for the Green's function G, i.e. by 
expanding G in powers of Gq. Yet, the corresponding 
expansion for the quasi-stationary state e nerg y converges 
very slowly. For that reason, we proposed a more effi- 
cient expansion scheme in which G is expanded in powers 
of the Green's function G^{E) = {E - K - W)-\ cor- 
responding to the outer potential W. From Eq. Q it 
immediately follows that 



G 



+ G^{U- Uo) G 



1 + WG 



(6) 



Iterating Eq. ((HJ in powers of G-^ and then substituting 
the result into ||2J), one finds the desirable perturbative 
expansion for the energy and the width of the resonance. 
By truncating this series, one obtains the following first- 
order relation valid for the isolated metastable state: 

E = Eo + {^o\W\<^o) + ($o|M^G^(S)M/|$o). (7) 

The above equation can be solved iteratively for E = 
Er — iT/ 2 hy assuming that the energy shift A = Er — Eq 
and the width F are small compared to Eq and Vq — Eq. 
In such a case, one can put G^^(i?) ~ Gy^^iEo), thus 
reducing Eq. (0 to 

E = Eq + {^o\W\<i>o) + {<!>q\WG^{Eq)W\^o)- (8) 

By using the Schrodinger equation for G^^ one finally 
obtains the TPA expressions 0, ^| for the width F and 
for the energy shift A ~ E^cs — Eq of the quasi-stationary 
state: 



mk 



2mk 



[aXkiR) + x'kiR)] [aXkiR) 



(9) 

x'kiR)], (10) 



where k = V2mEo, a = ^/2m{V{R) - Eq), Xkir) = 

5R(Xfc^ ''(»')), and Xk (xi^^) stands for the irregular (out- 
going) solution of the Schrodinger equation for the outer 
potential. 

It follows from Eqs. © and fTn|l that both F and A are 
given in terms of bound and scattering state wave func- 
tions. Thus TPA essentially simplifies the treatment of 
tunneling, because the standard approximation schemes 
can be used for evaluation of $0 and Xk ■ For instance, by 
applying the semi-classical approximation, one obtains 
the improved Gamow formula for F p , [Tol | , which is use- 
ful for different applications [T3. In particular, an 
extension of Eqs. © and 111)11 to the multi-dimensional 
case can be found in Ref. 
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III. CORRECTIONS TO TPA AND THE 
CHOICE OF THE SEPARATION RADIUS 

The accuracy of Eqs. © and (|10|l can be determined 
by evaluating the leading correction terms. There are 
two types of corrections to TPA: (a) those due to the 
replacement of G by Gy^ in Eq. leading to Eq. and 
(b) those due to the replacement of Gy^{E) by G^y(£'o) 
in Eq. Q leading to Eqs. Q, ((Tn|l . The correction terms 
of the first type, (Ar)i, can be obtained by iterating JB)). 
One finds from the first iteration [^: 

(Ar)i ^ mV'jR) 

r l&[2m{Vo-Eo)Y/^' ^ ' 

Equation might suggest that an optimal choice of 
the separation radius corresponds to V'{R) = 0, i.e., the 
maximum of V (r). However, it has been demonstrated 
numerically [T5l ll6| that if the top of the barrier is close 
to the closing potential, such a choice is not optimal, since 
in this case Eqs. © and H10() become less accurate. The 
reason is that the energy shift A becomes appreciable so 
that O cannot be replaced by ©. 

This can be illustrated by considering a square- well po- 
tential discussed in Ref. [ij: V(r) = l{l + l)/2mr'^ - 
for r < Ri, and V{r) — l{l ~\- l)/2mr^ for r > where 
the top of the barrier coincides with the closing potential 
(see Fig. 121). For a P-wave resonance, the numerical cal- 



V(r) 




FIG. 2: Metastable state in the square-well potential of 
Ref. [l7l |. The potential parameters (C/o=-51.6 MeV and 
i?i=6 fm) correspond to a P-wave resonance at E^cs ~ 1 
keV and F =55.9 eV. The bound state in the potential U{r) 
lies at iJo=275 keV. The top of the centrifugal barrier is at 
V(fli)=1.16 MeV. 

culation gives E'rcs = 1 keV and F = 55.9 eV ^3|- Now 
we apply the TPA by taking the separation radius at the 
boundary, R = Ri. The corresponding inner potential 



U{r) has a bound state at i?o=275 keV. However, the 
corresponding energy shift, A=— 300 keV, is of the same 
order of magnitude as the energy Eq of the bound state. 
Consequently, the replacement of G^{E) by G^y(£'o) in 
(UI) leads to large corrections to the resonance energy and 
the width so that Eqs. © and Hl()(l cannot be used. 

Fet us estimate the correction term (AF)2 = T{Eq + 
A) — F(-Eo) due to such a replacement. One can use the 
semi-classical Gamow formula, F oc exp(— 2 J^^ \p[r)\dr, 
with ri and r2 being the inner and outer classical turn- 
ing points, respectively, and \p{r)\ = y^2m[V{r) — Eq]. 
Approximating V{r) for ri < r < r2 by the inverted 
harmonic oscillator, one obtains 

(AF)2 . p V2^ , V2^(r2-ri)^ 
— - — ~ A / — , =ar sa tt _ -A , 

r Jr^ ^V{r) - Eo - Eo 

(12) 

where V= max V{r) and A is is given by Eq. (|10(l . 

Thus, in order to reduce the correction term (AF)2, 
one needs to minimize the energy shift A. It follows 
from 1)1 0|l that A can be strongly suppressed by taking 
the separation radius R deeply inside the barrier. Indeed, 
Eq. (|1U|) contains a product of regular and irregular wave 
functions, which do not vary considerably under the bar- 
rier. However, the factor |$o(i?)P decays exponentially 
with R. Therefore, by taking the separation radius R far 
away from the boundary, i? 3> ri, one finds that A — s- 0, 
and Eq — > E^cs- As a result, can be replaced by 0, 
leading to Eqs. © and (fTn|l . 

To illustrate this point, let us again consider the exam- 
ple of a P-wave resonance discussed above. By taking the 
separation radius R > Ri, we readily find that A — > 
as i? — i?i increases, and Eq Ej-^s— 1 keV. The result- 
ing value of F=55.3 eV is very close to the Gamow-state 
value of the resonance width. The separation radius R 
cannot be chosen too close to the outer classical turning 
point since in such a case V{R) — > Eq and the correction 
term 1)11(1 becomes important. In fact, Eqs. 1)11(1 and 1(12(1 
define the lower and upper limits of R. 

In the following, we discuss the TPA results for proton 
and neutron resonances in the realistic average nuclear 
potential. The approximate TPA expressions are com- 
pared to the resonant states results obtained using the 
GAMOW code [3. 

A. Comparison with Gamow-state calculations 

Consider single-nucleon resonances in the Woods- 
Saxon (WS) potential Vwsi'i') represented by a sum 
of central, spin-orbit, centrifugal, and Coulomb terms. 
Here, we apply the parametrization of Ref. jisf. namely: 
i?o=l-17A^/^ fm, a — 0.75 fm for the central term, and 
= 0.2C/o and = I.OIAVs for the spin-orbit po- 
tential. We calculate A and F according to Eqs. © and 
((10(1 by varying the separation radius R inside the bar- 
rier, starting with the barrier radius, R, corresponding 
to the maximum of V{r). 
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TABLE I: TPA calculations for the O/iii/2 and Isxji pro- 
ton resonances with energy i5ros=l-5MeV in a WS potential. 
The calculated widths Ftpa and the corresponding correc- 
tions mm and 11211 are shown relative to Ftpa for several 
values of the separation radius R. The actual accuracy of the 
TPA, Ar/r, is given in the last column. If the correction is 
marked zero, it means that it is below 0.1%. 



R-R 


A 


Ftpa 


(AF)i 


(AF)2 


(AF) 


(fm) 


(keV) 


(MeV) 


Ftpa 


Ftpa 


F 


Oftii/2 Gamow state: Ercs = 


=1.5 MeV, F =4.918 E-18 MeV 





-1.9 


4.931 E-18 





-1% 


-0.3% 


1.59 


-0.27 


4.919 E-18 





-0.15% 





4.28 


-5.0 E-3 


4.919 E-18 











2Si/2 


Gamow state: Ercs = 


1.5 MeV, F =6.695 E-14 MeV 





-3.3 


6.727 E-14 





-2.3% 


-0.5% 


3.48 


-0.11 


6.709 E-14 





-0.1% 


-0.2% 


8.05 


-1.2 E-3 


6.746 E-14 


-0.7% 





-0.7% 



TABLE IL Same as in Table HJ except for the neutron Oii3/2 
and 1/5/2 resonances. 



R-R 
(fm) 



A 

(keV) 



Ftpa 
(MeV) 



(AF)i 



Ftpa 



(AF)2 



Ftpa 



(AF) 



Oi 



13/2 



Gamow state: £^08=1 MeV, F =1.834 E-6 MeV 




3.45 
5.96 



-11.9 
-0.47 
-6.1 E-2 



1.869 E-6 
1.844 E-6 
1.847 E-6 





-1.2% 
-1.5% 



-3% 
-0.1% 




-1.9% 
-0.5% 
-0.7% 



1/5/2 Gamow state: £;i.es=l MeV, F =9.271 E-2 MeV 




2.05 
4.05 
4.57 



-109 
-51.4 
-7.2 
+2.7 



1.227 E-1 
1.089 E-1 
9.876 E-2 
9.408 E-2 





-10.1% 
-12% 
-17% 



-25% 
-12% 
-1.7% 
+0.6% 



-32% 
-17% 
-6.5% 
-1.5% 



We begin with the high-^ narrow proton resonance 
Ohii/2. The parameters of the WS potential are appro- 
priate for "'^'^^Tm, which is a proton emitting nucleus. 
The potential depth, [/o=— 61.8823 MeV, was adjusted 
to the energy i?rcs=l-5 MeV. The resulting barrier ra- 
dius is .R=8.54fm {V{R)^17.U MeV) and the inner and 
outer turning points are ri=6.33fm and 7'2=71.15fm, re- 
spectively. Since E^cs ^ ^{R), the calculated Ohii/2 res- 
onant state has a very small width, r=4.918 10~^^ MeV. 
The results of TPA calculations are shown in Tabic I 
for different values oi R > R, together with the cor- 
responding correction terms to TPA: (AT) 1 /Ftpa lfTT|l 
and (Ar)2/rTPA Table I also displays the actual 

accuracy of the TPA, (Ar)/r, where (AF) = P - Ptpa- 
Since R > ri, the energy shift A is small for R — R. 
Therefore, the results of TPA are in a good agreement 
with the resonant-state calculations already for R = R. 
Next we consider the \ow-£, broader 2si/2 resonance at 
i?i.cs=l-5MeV, which is considerably closer to the top of 
the barrier V{R) = 9.43 MeV (.R=9.34fm). As shown 
in Table I, also in this case Ptpa nicely agrees with the 
numerical result, and the accuracy of the TPA is well 
estimated by Eqs. ((TT|l and ((T^ . 



Table II displays the TPA results for the Oii3/2 
and 115/2 neutron resonances in ^^'^Sn at an energy 
£"103=1 MeV. Here P is much larger due to the absence 
of the Coulomb barrier. As in the proton case, there is 
very good agreement with numerical calculations, pro- 
vided that R is taken far away from the turning points, 
inside the window determined by (|ll|l and 1)12(1 . and the 
results of TPA weakly depend on the separation radius 
R. This suggests that the separation radius can be elim- 
inated altogether from the TPA expressions. As demon- 
strated in the following section one can indeed modify 
the TPA in such a way. 



IV. MODIFIED TWO-POTENTIAL APPROACH 



A tunneling potential can always be written as a sum of 
attractive and repulsive parts, V(r) = Vattir) + Vrepir), 
where Vrep {r) becomes dominant at distances beyond the 
barrier radius (see Fig.|31). Therefore, starting with some 




FIG. 3: For r > f the tunneling potential V{r) can be approx- 
imated by its repulsive part Vre.p{r). R and r2 denote the TPA 
separation radius and the outer turning point {V{r2) = Eq), 
respectively. 



radius f, the total potential V{r) can be well approxi- 
mated by the repulsive component only. The value of f 
should be chosen in such a way that the attractive (e.g., 
nuclear) part can disregarded with a desired accuracy 
77 < 1: 



\l-Vrep{r)/V{r)\<r^ 



for r > r . 



(13) 



For instance, in the case of a square-well potential of 
Fig. 12 Eq. 1(13(1 is satisfied for any 77 provided that f>Ri. 
In most cases, the attractive potential rapidly decreases 
beyond the barrier radius, so that f is closer to R than 
to the separation radius R in the TPA (cf. Fig.|3|). 

Usually, the repulsive part Vrep{r) is well known, as 
well as the two linearly independent (regular and ir- 
regular) solutions Fk{r) and Gk{r) of the corresponding 
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Schrodinger equation. For instance, if Vrep{r) is a sum 
of Coulomb and centrifugal potentials, then -F'fc(r) and 
Gk{r) are the standard Coulomb functions. This implies 
that any solution of the Schrodinger equation with the 
potential V{r) can be written for r > f as the linear 
combination of Fk (r) and Gk (r) . 

Consider the bound-state wave function ^o{r) of the 
inner potential U{r) of Eq. Q. Since U{r) ~ Vrep{r) for 
f < r < R, ^a{r) can be expanded in this region as 



Mr) ^ ciGkir) + C2Fk{r) , 



(14) 



where k = (2rn£'o)^/^. The coefficients ci.2 and the en- 
ergy Eq are obtained from matching of the logarithmic 
derivatives at r = f and r = R: 



c^Gm+C2Fl{R) 
ciGk{R)+C2Fk{R) 
ciG",(f)+C2F^(r-) 



= -\p{R)l 



ciGkif)+C2Fk{f) $o(0' 



(15a) 
(15b) 



Note that $o(r) « exp[-|p(i?)|(r - i?)], for r > R. Solv- 
ing Eqs. (|15|l one easily finds 



Mr) Gkif) 



I + aa^' - f2{l + dk') 



(16) 



where a = \p{R)\, Uk = Fl^{R) / Fk{R) , 
G'^{R)/Gk{R), and 



ak 



h 



G'k{R)F'kir) 
F'kiR)G'k{r) 



f2 



G'k{R)Fk{r) 
niR)Gk{f) 



(17) 



The wave functions Fk{r) and Gk{r) are of the same 
order of magnitude in the asymptotic region, r ^ r2. 
However, in the classically forbidden region the regular 
wave function exponentially decreases and the irregu- 
lar one exponentially increases with decreasing r. Us- 
ing the semi-classical approximation, one can estimate 
Gk{R)Fk{f) ~ exp[— /i* |j)(r)|dr]. Therefore, the coef- 
ficients /i^2 are of the order of exp[— a(i? — f)], so the 

corresponding terms /i,2(l + a[. ^^) in Eq- (|1B|) are expo- 
nentially suppressed and can be neglected. As a result, 
the matching condition ll()ll can be written as 



%{f) _ G'^{f) 
Mr) Gk{f) ■ 



(18) 



The above equation constitutes the MPTA condition 
for the resonance energy E-^cs — Eq. In contrast to 
Eq. Hlt)|) . the relation H18|) does not exhibit any explicit 
i?-dependence. Consequently, there is no need to evalu- 
ate the bound-state wave function at large values of R 
well inside the barrier. Note also that for a narrow res- 
onance the irregular wave function Gk{r) is proportional 
to the real part of the outgoing (Gamow) solution: 



^out 



(r) cx Gfc„, (r) + zFfe^„ (r) = Ok^^^ (r), (19) 



with complex k^es = k — and fc^g^ = 2m{Eres — *r/2). 
Since for small T and 7 the imaginary parts of kres and 
Gk^cB (r) can be neglected, kres ~ k and Gk^^^ {r) ~ Gk{r) 
(also for the outgoing Coulomb wave function Ofc^ss {r) = 
Ofe(r)). In this case Eq. H18() represents the matching 
condition for the inner (bound state) wave function with 
the real part of the non-normalized outgoing wave. If 
the imaginary parts are not negligible, Eq. (|18|l can be 
still interpreted in terms of the standing-wave boundary 
condition at f, which means that the scattering phase 
shift is 7r/2. The requirement that the phase shift is 
equal to 7r/2 represents an alternative definition for the 
position of a resonance in the absence of a non-resonant 
phase shift. 

Consider now Eq. for the width. Since R > r, the 
outer wave function Xk{r) in the region r > R can be 
represented by the linear combination of the regular and 
irregular solutions of Vrep (r) , and the corresponding coef- 
ficients in the linear combination of Fk{r) and Gk{r) are 
directly related to the (non-resonant) scattering phase 
shift for the outer potential W{r). One easily finds 

Xk (r) = cos 4 Fk (r) + sin 4 Gk (r) for r > i? , (20) 

where the phase shift Sk is obtained from matching of 
logarithmic derivatives at the separation radius R: 



tan 6k 



Fk (R) f Ok - a 



GkiR) \ak 



a 



(21) 



Here we neglected the terms ~ exp(— 2ai?). Substituting 
(|14|l and (|20|l into ^ and taking into account the Wron- 
skian relation between Fk{r) and Gk{r) one obtains: 



r = cos 6k — 
m 



r*o(r)i 


2 r 


[Gk{f)\ 





(ofc 4- a){ak - a) 
(ofe -I- a){ak - a) 



(22) 



Note that in the classically forbidden region 
Fl,{r)/Fk{r) « \p{r)\ and G',{r)/Gk{r) « -\p{r)l 
so that the second term in brackets of H22|) can be 
neglected. In addition, cos 4 ~ 1, as follows from 
(|21|l . As a result, one arrives at the following simple 
expression for the width: 



Mr) 



Gk{f) 



(23) 



Thus, similar to Eq. (fTH|l for the resonance's energy, the 
separation radius R does not appear explicitly in the ex- 
pression for the width. 

Equations (|18|) and H23|l represent the final result of the 
modified two-potential approach (MTPA). Despite their 
simple appearance, these expressions are very accurate. 
In fact, the accuracy of the MTPA is practically the same 
as that of the TPA since the former was derived from the 
latter by neglecting only small correction terms of the 
order of the accuracy of the TPA itself. For instance, for 
the previously discussed case of the P-wave resonance 
in the square well potential, one finds £'ics=lkeV and 
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r=55.3eV, i.e. the same result as in TPA. In general, 
the accuracy of the MTPA can be estimated by means 
of the parameter ry, which defines the lower limit for the 
matching radius f . However, one has to keep in mind that 
f cannot be very large since the derivation of Eqs. H18() 
and 1)23(1 is valid only for a{R — r) ^ 1. Therefore, the 
value of T] is restricted by Eq. in which R is replaced 
by f. 

It is worth noting that an expression similar to Eq. H23() 
was used in Refs. 0,113 for calculating partial widths for 
proton emission. The corresponding formula that applies 
to the single channel case can be written as 



r(r) = 



3? 



Okir) 



(24) 



where r is large. The r-dependence of r(r) is weak but 
it can be reduced if one takes a more appropriate expres- 
sion: 



r(r) = A 

m 



\YjOut 



(r) 



(25) 



Another expression for the width can be derived from the 
continuity relation for the resonant states 7, 21]: 



r(r) = 



I -q/ouu* ^^-j.q,out (|.) _ vi/o«t'(7-)\i/o«t*(,.) 
2m 



(26) 



This form is completely independent of r in a wider range 
pH |2^ and furnishes a value which is equal to that com- 
ing from the imaginary part of the energy. However, for 
very narrow resonances, it is difficult to calculate the 
imaginary part of the energy with sufficient precision. 
The expression (I23II derived in the MTPA replaces the 
Gamow wave function with the (real) bound-state wave 
function $o- Finally, let us emphasize that while Ea. H24|) 
resembles the MTPA expression, it is based on different 
approximations and boundary conditions. On the other 
hand there are close connections between the MTPA and 
the R-matrix theory, also employing real-energy eigen- 
states, see Sec. Ivl 



A. MTPA: numerical examples 

We present below in Table HI the results of the MTPA 
for the widths H23|) of resonances discussed in Tables I 
and II in the context of TPA. (Since the MTPA reso- 
nance energies are very close to the exact result, they are 
not displayed.) One finds that the MTPA reproduces 
the width almost with the same accuracy as the TPA, 
provided that the matching radius r is large enough to 
ensure that the contribution from the nuclear attractive 
potential is small (ry <C 1). It follows from Table HI that 
ij controls the accuracy of MTPA rather well, except for a 
broad neutron resonance I/5/2 when AT /T reaches 10% 
at f = R + 2.96. In this case, the matching radius is 



TABLE III: Similar as in Tables Q and ^ except for the 
MTPA. 



f-R 
(fm) 


V 


Pmtpa 
(MeV) 


(Ar) 
r 


O/111/2 Gamow state: Bros =1.5 MeV, T =4.918 E-18 MeV 


0.17 
1.55 
3.09 


0.11 
0.02 
0.003 


4.665 E-18 
4.87 E-18 
4.909 E-18 


5% 
1% 
0.2% 


2si/2 Gamow state: _E,.os=1.5 MeV, V =6.695 E-14 MeV 


0.26 
1.66 
3.24 


0.06 
0.01 
0.001 


6.577 E-14 
6.675 E-14 
6.692 E-14 


2% 
0.3% 



0*13/2 Gamow state: _B,.os = lMeV, T =1.834 E-6 MeV 


0.18 
1.45 
2.99 


0.15 
.04 
.007 


1.736 E-6 
1.814 E-6 
1.831 E-6 


5% 
1% 
0.1% 


1/5/2 Gamow state: £'rcs=lMeV, T =9.271 E-2 MeV 


0.18 
1.38 
2.96 


0.13 
0.035 
0.005 


8.998 E-2 
8.856 E-2 
8.373 E-2 


3% 
4% 
10% 



quite far away from the barrier radius, so that the accu- 
racy of the MTPA is given by Eq. ((TT]) (with R replaced 
by f). This is well confirmed by Tabic II, which shows 
the corresponding correction term. 



V. COMPARISON WITH THE R-MATRIX 
THEORY 

It is interesting to compare the final expressions of the 
MTPA, Eqs. JTHl and with the results of the R- 

matrix theory In the latter method, the space is di- 
vided into internal and external regions by a hard sphere 
of the radius ac, and a complete set of the internal wave 
functions u\(r) is introduced, J^" u\{r)u\'{r)dr = d\x'- 
The internal wave functions obey real boundary condi- 
tions for the logarithmic derivative: 



u\{ac) 



B. 



(27) 



The value of B determines the R-matrix energy eigenval- 
ues E\. It is convenient to choose B so that one of the 
eigenvalues would coincide with the position of the reso- 
nance where the value of the phase shift is equal to 7r/2. 
This defines the so-called "natural" boundary condition: 



<(Qc) 



(28) 



In this case, one can apply the one-level approxima- 
tion, i.e., approximate the resonance with a single Breit- 
Wigner term. The corresponding width Pa becomes: 



(29) 
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where g\ denotes the energy derivative of the level shift 

One finds that Eqs. 10 and ^ of the MTPA for- 
mally resemble Eqs. (I28II and (|29|l by choosing Qc — r 
and taking g\{ac) — 0. Yet, the inner wave function 
^i){r) of the MTPA is different from the internal wave 
function u\{r) of the R-matrix theory. The latter is to- 
tally confined inside the inner region, whereas '^o{r) is a 
true "bound state" wave function of the inner potential 
U{r). Therefore, their normalizations are different. 

The essential problem of the R-matrix theory is a 
proper choice of the matching radius Cc, which remains 
a free parameter. A different choice of Oc does affect 
the results of R-function calculations in a one-level ap- 
proximation. Moreover, there should exist an optimal 
matching radius, for which the results of the R-matrix 
calculations are close to the exact results However, 
except for some simple cases (e.g., the square well po- 
tential), the optimal matching radius cannot be simply 
prescribed. In contrast, the MTPA is not sensitive to 
the matching radius f, provided that it is taken inside 
the "window" defined by Eqs. ifTTjl and (fT^ . This is an 
essential advantage of MTPA over the R-matrix method. 
(For critical discussion of the R-matrix expression for the 
resonance width, see Refs. [8Ll2lLl2^.) 

VI. SUMMARY 

This paper contains a detailed investigation of the two- 
potential approach to the one-body tunneling problem. 
It has been found that TPA becomes extremely accurate 
if the separation radius, dividing the entire space into 
the inner and the outer regions, is taken deeply inside 
the barrier, but not too close to the outer classical turn- 



ing point. From a minimization of the leading correction 
terms, we obtained simple expressions for the upper and 
lower limits of the TPA separation radius. The high ac- 
curacy of the method was demonstrated explicitly by a 
detailed comparison with Gamow resonances of a realis- 
tic nuclear potential. 

Furthermore, we have found that the TPA can be fur- 
ther simplified by taking into account the properties of 
regular and irregular solutions of the Schrodinger equa- 
tion under the barrier. The final expressions of the mod- 
ified two-potential approach formally resemble those of 
the R-matrix theory with the "natural" boundary condi- 
tions. However, the internal wave function of the MTPA 
is considerably different. In addition, contrary to the R- 
matrix theory, the corresponding matching radius of the 
MTPA is well defined. This makes MTPA particularly 
suitable for practical applications. 



VII. ACKNOWLEDGMENTS 

One of the authors (S.G.) is very grateful to E. Vogt 
for inspiring discussions and important suggestions. S.G. 
also thanks TRIUMF, Vancouver and JIHIR for kind 
hospitality. This work was supported in part by the 
U.S. Department of Energy under Contract Nos. DE- 
FG02-96ER40963 (University of Tennessee), DE-AC05- 
00OR22725 with UT-Battelle, LLC (Oak Ridge National 
Laboratory), and DE-FG05-87ER40361 (Joint Institute 
for Heavy Ion Research). T.V. acknowledges support 
from the Hungarian OTKA Grants No. T029003 and 
T037991. The support from the University Radioactive 
Ion Beam consortium (UNIRIB) is gratefully acknowl- 
edged. 



[1] F. Hund, Z. Physik 43, 805 (1927). 

[2] L. Nordheim, Z. Physik 46, 833 (1927). 

[3] J. R. Oppenheimer, Phys. Rev. 31, 66 (1928). 

[4] G. Gamow, Z. Physik 51, 204 (1928); 52, 510 (1928). 

[5] R.W. Gurney and E.U. Condon, Nature 122, 439 (1928); 

Phys. Rev. 33, 127 (1929). 
[6] A.M. Lane and R.G. Thomas, Rev. Mod. Phys. 30, 257 

(1958). 

[7] J. Humblet and L. Rosenfeld, Nucl. Phys. 26, 529 (1961). 
[8] A. Arima and S. Yoshida, Nucl. Phys. A219, 475 (1974). 
[9] S.A. Gurvitz and G. Kalbermann. Phys. Rev. Lett. 59, 
262 (1987). 

[10] S.A. Gurvitz, Phys. Rev. A 38, 1747 (1988). 

[11] S.A. Gurvitz, Michael Marinov Memorial Volume, Mul- 
tiple facets of quantization and supersymmetry, p. 91, 
Eds. M. Olshanetsky and A. Vainshtein (World Scien- 
tific, 2002). 

[12] D.F. Jackson and M. Rhoades-Brown, Ann. of Phys. 105, 
151 (1977). 

[13] T. Berggren and P. Glanders, Nucl. Phys. A473, 189 
(1987). 



[14] B. Buck, A.C. Merchant, and S.M. Perez, Phys. Rev. 
Lett. 65, 2975 (1990); B. Buck, J.C. Johnston, A.C. Mer- 
chant, and S.M. Perez, Phys. Rev. C53, 2841 (1996). 

[15] S. Aberg, P.B. Semmes, and W. Nazarewicz, Phys. Rev. 
C56, 1762 (1997). 

[16] P. Talou, N. Carjan, C. Negrevergne, and D. Strottman, 
Phys. Rev. C62, 14609 (2000). 

[17] E. Vogt, Phys. Lett. B 389, 637 (1996). 

[18] T. Vertse, K.F. Pal, and Z. Balogh, Comp. Phys. Comm. 
27, 309 (1982). 

[19] E. Maglione, L.S. Ferreira, and R.J. Liotta, Phys. Rev. 

Lett. 81, 538 (1998); Phys. Rev. C59, R589 (1999). 
[20] E. Maglione and L.S. Ferreira, Phys. Rev. C 61, 047307 

(2000). 

[21] B. Barmore, A.T. Kruppa, W. Nazarewicz, and T. 
Vertse, Phys. Rev. C62 054315 (2000). 

[22] T. Vertse, A.T. Kruppa, B. Barmore, W. Nazarewicz, 
L. Gr. Ixaru, and M. Rizea, in Proceedings of the In- 
ternational Symposium on Proton-Emitting Nuclei, Oak 
Ridge, edited by J. Batchelder, AIP Conf. Proc. No. 518 
(AIP, New York, 2000), p. 184. 



8 



[23] G. Michaud, L. Scherk, and E. Vogt, Phys. Rev. C 1, 864 Emitting Nuclei, Oak Ridge, edited by J. Batchelder, AIP 

(1970). Conf. Proc. No. 518 (AIP, New York, 2000), p. 173. 

[24] A.T. Kruppa, P. Semmes, and W. Nazarowicz, in Pro- 
ceedings of the International Symposium on Proton- 



